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Abstract 

We provide a semiclassical description of the electronic transport through graphene n-p junctions 
in the quantum Hall regime. A semiclassical approximation for the conductance is derived in terms 
of the various snake-like trajectories at the interface of the junction. For a symmetric (ambipolar) 
configuration, the general result can be recovered by means of a simple scattering approach, pro- 
viding a very transparent qualitative description of the problem under study. Consequences of our 
findings for the understanding of recent experiments are discussed. 



PACS numbers: 73.22.Pr, 73.43.Jn, 03.65.Sq, 73.23.Ad 



Graphene is a two dimensional carbon based material which, due to its remarkable elec- 
tronic and mechanical properties as well as potential applications, has been the subject of 
intense research activity in physics, chemistry and material sciences [H [2] . At the origin of 
this interest is the observation that the electronic low energy dynamics in graphene is gov- 
erned by a Hamiltonian very similar to that of the 2-dimensional relativistic Dirac equation. 
This has a number of remarkable consequences, two of the most striking being an anomalous 
quantum Hall effect [H E] and the existence of Klein tunneling jSHTj. 

Experiments where both quantum Hall physics in graphene and Klein tunneling are 
at play were pioneered by Williams and collaborators, who measured the conductance of 
graphene n-p junctions in a high magnetic field in the quantum Hall regime [Sj- This was 
possible due to the deposition of a metallic topgate partially covering the graphene sheet. 
By independently varying the applied top and backgate voltages the resulting electrostatic 
potential creates negatively and positively doped regions in the graphene sample. In this 
way n-n, p-p, and n-p junctions were produced. The two former did not show any striking 
feature. In contrast n-p junctions, where the transition from the electron region to the hole 
region of the sample has to take place through Klein tunneling, showed a quite unexpected 
behavior. Conductance plateaus at values Go/2 and 36*0/4 (Go = ^e^/h) were observed, 
at odds with the sequence {2n + l)G'o, n G Z, expected for the quantum Hall effect in 
graphene. Other experimental groups have observed these and other interesting plateaus in 
more elaborated set-ups P4TT]. 

These observations were explained by Abanin and Levitov [12] through a "quantum 
chaos hypothesis", which can be summarized as follows. At the interface between the n 
and p regions of the junction, the electrons experience a succession of Klein tunneling and 
skipping-orbit like propagation. It was suggested that the mode mixing caused by this 
mechanism possibly leads the probability to be transmitted or reflected in a given mode to 
be perfectly "democratic" . The Landauer-Biittiker formula for the conductance under this 
hypothesis gives G = G^NnNp/ {Nn + A'p), where Nn and Np are the number of channels in 
the n and p regions. This simple formula agrees with the observed conductance plateaus. 

As already noted [12j, this interpretation fails to provide a complete explanation of the 
experimental findings. Indeed, the "quantum chaos hypothesis" corresponds to the assump- 
tion that the scattering matrix describing mode mixing along the junction can be statistically 
modeled by random matrix theory. This hypothesis implies that the average scattering ma- 
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trix coefficients are equal. However, it also predicts universal-conductance-like fluctuations 
[131 [E], which are quite robust and are not suppressed by adding disorder [TSl US]- Ex- 
periments show reasonably well defined plateaus, but no significant fiuctuation amplitudes 
around a mean value. An alternative approach p!7l [18] showed that the edge boundary con- 
ditions affect the valley polarization of the zero energy Landau level. Clean edges produce 
new conductance plateaus in the ballistic regime, however not at the observed values. The 
graphene n — p junction plateaus remain an experimental mystery. The purpose of this 
Letter is to present a semiclassical analysis that provides a better understanding of mode 
mixing processes at the graphene n-p interface regions, contributing to dispel the puzzle. 

Mode mixing depends on several system-specific features: presence or absence of disor- 
der, steepness of the potential barrier, possible diffractive effects when the electrons transit 
from their skipping orbit motion along the sample edges to the snake-like motion along the 
junction interface, etc. To specifically focus on the mixing occurring at the interface, we 
consider the model sketched in Fig. [1} a perfectly clean graphene sample, the edges of which 
match smoothly with a step-like potential interface. The Landauer-Biittiker conductance G 
across the junction is expressed as 



where the sum runs over all the incoming skipping modes k, a {a is the valley index) in for 
instance the n region, and the T^^a are the corresponding transmission probabilities across 
the junction, for which we present a semiclassical description. 

It is useful to begin by analyzing what are the results expected from a simple classical 
point of view. Let us call > and Vp < the electrostatic potential in the n and p 
regions respectively (we assume the chemical potential ^ = 0), the applied magnetic 
field, A the corresponding vector potential, and vp the Fermi velocity. In the interior of 
either regions, electrons experience a cyclotronic motion of opposite direction and of radius 
Rn,p = \Vn,p\/{€.B;^VF). At the sample edge, the classical trajectories follow the skipping 
motion illustrated on Fig. [T] and can be characterized by the angle at which they hit the 
boundary. We call this angle On for trajectories in the n region, and it — 6p for trajectories in 
the p region (as it turns out to be more convenient to label the direction of 11 = p -|- eA(r), 
which is antiparallel to the velocity in the hole region). Between two successive bounces, 
trajectories progress a distance L„p = 2Rn^pSm9n,p along the edge. The motion is similar 
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FIG. 1: (Color online) Scattering along a boundary intercalated with a step-like n-p junction. The 
different symbols are explained in the text. 



at the n-p interface, except that at each encounter with the potential barrier there is a 
probability 

^ _ sin 9n sin 9p . 
~cos2(^) 

to be transmitted from the electron side to the hole side (or reciprocally) between two 
trajectories fulfilling the Snell-Descartes relation 

Rn cos 9n = Rp cos 9p . (3) 

Let us consider an incoming mode (k, a) from the n region. Semiclassically, the mode is 
built on a one-parameter family of skipping trajectories corresponding to some fixed 6'^'". 
The different trajectories within the family are labeled by the abscissa xq G [— L„(^^'"), 0] 
at which they last bounce before entering the interface region. The evolution of this set 
of trajectories under the above classical dynamics can be characterized by the probabilities 
Un{x) and Up{x) to emerge at x on the electron side (with angle ^^'°) or on the hole side 
(with the corresponding angle given by Eq. ([s])). If the width L of the junction is signifi- 
cantly larger than L„ and Lp, it can be shown that Un{x) and Up{x) converge both toward 
1/2 exponentially quickly with the number of bounces on the junction. As a consequence, 
since whether a trajectory is transmitted or reflected only depends on the side it emerges 
from after the last scattering at the interface, the classical probability of transmission is, 
for a large junction, given by the ratio Lp/(L„ + Lp). It should be stressed that this clas- 
sical transmission probability, even if averaged on the angle 9^'°' specifying the incoming 
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mode, does not correspond to what a "quantum chaos hypothesis" would require, namely 
a probability given by the proportion of available classical phase space on each side of the 
junction. 

Let us now turn to the description of the quantum transmission. For the sake of defini- 
tiveness we discuss the case of infinite mass boundary at the edge of the graphene sample. 
We expect no qualitative differences for either armchair or zigzag boundary conditions. Let 
us introduce the actions 



Semiclassically, the edge modes in the n region are built, within the leads, on the skipping 
trajectories bouncing on the boundary with an angle ^ fulfilling the quantization condition 



(again, a = ±1 is the valley index). Note there is no k = level for a = +1. 

We start our semiclassical discussion with the particularly simple case of ambipolar junc- 
tions, namely junctions such that Vp = —Vn- In that situation, the Snell relation (|3]) tells 
us that Qn = Op for any pair of angles, and therefore, Ln{6n) = Lp{6p). As a consequence, 
the different trajectories within the family {On = 6'^'",xo G [— 0]} follow a com- 
pletely independent history. The propagation of the corresponding amplitudes along the 
junction can be performed for each of them separately and is obtained from the successive 
multiplications of two 2x2 unitary matrices. The first one. 



with here On = Op = describes the propagation on the electron and hole side between 
two interactions with the interface. The ±7r/2 are the Maslov phases originating from the 
focal point met on these pieces of trajectory, and ^ = —On is a Berry phase. The second 
matrix describes the transmission and reflection between the electron and hole sides. It can 
be expressed as (in a form not restricted to the ambipolar case) 




(4) 




5„(^;;'") = 27ih{K - a/4) 



(6) 
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Length of the interface Length of the interface Length of the interface 

FIG. 2: (Color online) Left: Conductance (in units of Gq) as a function of the length of the interface 
L (in units of the magnetic length Ib = y^hJ{eBz)) for the case of 13 edge channels (filling factor 
7) in the ambipolar (thick black line) and generic (red dashed line) cases. Evidently, no plateau 
in the conductance is reached in either cases. Middle and right: Transmission of edge channels 
number 4 and 13, for which the local transmission probability in Eq. ^ is respectively 71 = 
0.87 and Tls = 0.19, for semiclassical (thick black line) and classical (red dotted line) dynamics. 
Substantial difference between both behaviors indicates strong influence of interferences between 
classical trajectories. 

with r„ = rp = -cos((6'„ + 6'p)/2)/cos((6'„ - 0p)/2) and tn,p = -ism9n,p/ cos((6'„ - 9p)/2). 
Note the phases e^^^"--^'^^"'^'^^'^ can be interpreted also as Berry phases. The transmission 
probability for a given trajectory is then given hj T{9^^'^,Xo) = 1(0, 1)(PD)^(1, 0)^1 with 

the integer part of {L + \xo\)/Ln, which therefore depends on xq. The total transmis- 
sion for the mode (/t, a) is then obtained as T^^a = J^j^ dx^TiO'^'"' ,xq). The resulting 
conductance after summation on the modes (k, a) is shown on Fig. |2} 

In the generic case of arbitrary Vn and Vp, a semiclassical expression for the conductance 
can be derived from the Fisher-Lee-Baranger-Stone equations [121 120] relating the transmis- 
sion and reflection coefficients T^^a and R^^a to the Green function Q{r" , r') and the incoming 
mode Xii,a{^')- our case the transmission can be written as 

T^,c, = I dy"tiAr")vFaJ^,air") , (9) 

Jy">0 

with 

tK,a(r") = -ih / dy'g{r",r')vF(T^XK,ai^') ■ (10) 

Jy'<0 

The same results holds for i?K,a except that the integral on y" should be taken for negative 
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values of y" (i.e. in the electron region). The incoming mode is semiclassically expressed as 

X.Jr) = 4=f^ E e^^^^-(^)-^ ..^^ ^ (11) 
VI sin ^(y) I ,t^i \^ e*2«(f) y 

with ^(y) the angle (> 0) such that y = i?„(cos^^'" - cose{y)), = {eBQRl/2){e{y) - 

sin(2^(?/))/2), and the normalization coefficient A^^a = (4f^i?„sin 

function we take the semiclassical approximation derived in [21] for the graphene Hamil- 
tonian, properly modified to account for Klein tunneling at the junction interface. The 
semiclassical Green function is expressed as a sum over all trajectories j joining r' to r" 

(12) 

where Sj is the action integral along the orbit j, Jj the stability determinant, fij the Maslov 
index counting the focal points met by the trajectory (counted negatively on the hole side), 
C,j is a Berry phase equal to half of the angle of rotation of the vector 11 = p + eA(r) 
(which is parallel to the velocity on the electron side and antiparallel to it on the hole side), 
= -^{Jie), V~ = ;^(^_i )) with 9 the direction of 11, and e(r) = +1 or —1 depending 
on whether r is on the electron or hole side. The only difference introduced by the presence 
of the potential barrier [22] is that bounces on the interface, indexed by 'jn,p, as well as 
transmission across it, indexed by rjn^p, are taken into account through the coefficients rn,p, 
tn,p introduced in Eq. ^ (in the definition of which the angles ^(T"'P'^"'p) are used). 

A typical skipping orbit fulfilling the quantization condition Eq. (|6| is displayed on 
Fig. [l} The orbit is characterized by the angle 9" at which it arrives at the Poincare 
section x = L, the number of excursions rUp on the hole side, and the number of traversals 
k of the interface. The number of excursions on the electron side is then mn{'mp,9") = 
[{W — L" — 2mpRp sin 9p) / (2_R„ sin 6'„)] with [■] the integer part, and L" = Rp{sm 9p — sin 9") 
for transmission and i?„(sin^„ + sin^") for reflection. The trajectory should also be char- 
acterized by the ordering of the various excursions in the n and p sides of the junction. 
However these various orderings will correspond to the same amplitude, and therefore just 
contribute degeneracy factor given by n{mp,k = 2k') = (Xi') (""^T"^"^') reflection. 



and Q{mp, k = 2k" + 1) = (™,^) (™"^^,T''^) for transmission. Performing the integral in Eq. (10) 
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in the stationary phase approximation, we get 



with 



T 

■J- K 



Rn 



0T 



2R n sin 9n J —0rp 

1 r<^R 



sm &n 



do" cose"n 



(13) 
(14) 
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mr,=0 



k" 



t^tp 



fi(mp,2fc" + l)(15) 



Mp(0") 



mr,=0 



ttV 

rnTpJ 



(16) 



where Mp{e") = [{W - L")/(2i?pSin0p)] and Or = min(^p, vr - Op), Or = min(e^'°, vr - ^;;'°). 
The corresponding curves for the conductance in Eq. ([T| are shown on Fig. [2] 

A few comments are in order. First, we stress that in the ambipolar case Vp = —Vn, the 
results obtained within the Baranger-Stone framework are, as expected, strictly identical to 
those derived earlier with the "scattering matrix" approach. Second, the equilibration and 
thus the saturation of the conductance expected at the classical level has to be contrasted 
with the persisting large oscillations of the semiclassical conductance. Finally and more 
surprisingly, we observe that the mean of the semiclassical prediction can differ significantly 
from the classical limiting value. This is illustrated, for instance, at the rightmost panel of 
Fig.|2j 

Understanding the oscillating pattern in the conductance as a function of the length of the 
interface in Fig. [2] is relatively straightforward in the ambipolar case, since the "scattering 
matrix" can be interpreted as that of a rotation on the Bloch sphere, acting on the vector 
composed of the reflection and transmission probability amplitudes of a given mode. In 
particular convergence to a given transmission probability, corresponding to a vector pointing 
at a fixed latitude, cannot be achieved. Outside of the ambipolar case, this interpretation 
is no longer valid but the pattern of the conductance in Fig. [2] is similar enough to let us 
believe this picture remains qualitatively correct. 

The advantage of a semiclassical description, compared for instance to an exact numerical 
calculation of the conductance using recursive Green function techniques, is that it makes 
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it possible to discuss the expected consequence of various modifications of the model we 
consider. For instance we do not expect that including a finite width in the potential 
step (as long as <^ Ib) or changing the edge boundary conditions will qualitatively 
modify the oscillating pattern of the conductance. In the first case, the local probability 
transmission Eq. ^ will be somewhat smaller, actually improving the speed of classical 
convergence. In the second case, the quantization condition in Eq. (|6| will be modified [23] , 
but without affecting the basic mechanism at play here. A geometry closer to the one used 
in experiments, i.e. with the junction perpendicular to the edge of the ribbon, would imply 
some diffraction at the edge-junction corner (this aspect will be discussed in more details in 
|24j). This, as more generally the inclusion of a weak disorder, would somewhat diminish the 
amplitude of the conductance oscillations and bring the average of the semiclassical results 
closer to the classical conductance. 

As increasing the amount of disorder could only bring the system toward a chaotic limit, 
characterized by universal conductance fluctuations around a mean given by the classical 
(democratic) expectation, our study indicates that within a model of perfectly coherent 
electrons, the intrinsic properties of the n-p junction cannot produce the experimentally 
observed plateaus. These considerations suggest the existence of inelastic processes occurring 
in the vicinity of the junction, possibly reducing the coherence length to values smaller 
than L. Further experiments, varying the ratio l^/L, are expected to provide further insight 
on this issue. 

We acknowledge a fruitful discussion with Alfredo Ozorio de Almeida. This research was 
supported by the CAPES/COFECUB (project Ph 606/08). 
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